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Abstract 

We briefly review the properties of radially growing interfaces and their connection to biological 
growth. We focus on simplified models which result from the abstraction of only considering 
domain growth and not the interface curvature. Linear equations can be exactly solved and the 
phenomenology of growth can be inferred from the explicit solutions. Nonlinear equations pose 
interesting open questions that are summarized herein. 
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I. INTRODUCTION 



Since the origins of the interdisciplinary exchange among mathematics and biology, the 
development of the form and structure of living organisms has always been considered a 
fascinating topic l|, [2J. Already present in these seminal works, the necessity of considering 
physical and mathematical laws, together with the relevant biological principles, has been 
increasingly recognized in the scientific literature along the years. One such example is the 
examination of the properties relating to the architecture of cell colonies, as already noted 
by Murray Eden jjj, 4]. The overall appearance of a living organism is conditioned by its 
genetic expression. Indeed, it escapes to nobody that twins are similar to each other as whole 
individuals. However, when we regard smaller scale properties of them, like the dermal ridge 
count of fingerprints or the patterns of retinal venation, their differences are greater even in 
the case of monozygotic twins. Differences that are expected to grow as our reference spatial 
scale decreases, and that are particularly evident in the architecture of cell colonies. These 
and other observations unveil that morphogenetic processes are not completely determined 
by the genetic expression. Of course, it is widely recognized nowadays that environmental 
factors play a determinant role in morphogenesis. Both aside and as a consequence of this 
one would be interested in determining how random factors affect the process of growth, 
and so how is this process constraint by the laws of probability. 

The Eden model was introduced in order to shed some light into this question 0, 4|. It 
is a simplified model for biological growth, which results from the abstraction of neglecting 
many of the real aspects of cell colonies development. It concentrates just in the appearance 
of new cells in the colony periphery; once introduced, cells are never removed from it. This 
model is to be considered in (Z^)^ ' 1 ^, where d EN denotes the spatial dimension, and and 
1 stand for an empty and an occupied site respectively. Starting from a single cell at the 
origin, and following a set of probabilistic rules which dictates the frequency and manner in 
which new cells are introduced in the colony [si, Q], a radial form develops in the long time. 
Its macroscopic shape is affected by the underlying lattice structure, and its interface shows 
fractal properties that are independent of it 5, (|. 

Although a rigorous connection have never been established, it is widely accepted that 
the fluctuations of the Eden interface can be described with some suitable stochastic partial 
differential equation 

Baa 

. Such a relation has been proposed in light of numerical results 
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derived with the so called cylindrical Eden model Q. In this case, for d = 2, the system is 
a strip of infinite length and finite width L, and the initial condition is a whole semistrip. 
For large enough L and periodic boundary conditions the system obeys the Family- Vicsek 
i. e., the two-points correlation function is of the form 



ansatz 

(h(x, t)h(x', t)) = \x- A 2a f (^JT 1 ) . (!) 

where (•) is the average over a large number of realizations, /(■) is the scaling function, and 
a and z are universal quantities known as critical exponents: a is known as the roughness 
exponent and z as the dynamical exponent, which measures the velocity at which the cor- 
relations travel. The function h(x,t) expresses the height of the interface with respect to 
the initial condition at some given position x and time t. The ratio of the critical exponents 
(3 = a/z constitutes a new exponent measuring the velocity at which the interface width 
increases during the first stages of growth . For the dimensionality considered the sim- 
ulations have measured a « 1/2 and z ~ 3/2, which place this discrete model in the same 
universality class as the continuum equation 

d t h = vV 2 h+-{Vhf + F + £(x,t), (2) 



which is known as the Kardar-Parisi-Zhang (KPZ) equation [9|. In this case £(x,t) is a 
Gaussian distributed spatiotemporal noise delta correlated in both space and time, and u, 
A and F are real positive parameters. Of course, the cylindrical Eden model differs from its 
radial counterpart in two features: the interface is curved and grows laterally. Herein we will 
consider an abstraction of this problem and focus exclusively on lateral growth. As discrete 
models are usually placed in the universality classes defined by continuum equations, we will 
study the dynamics of such equations defined on uniformly growing domains. 



II. LINEAR GROWTH 



In this section we will summarize some of the recent results obtained with linear equa- 
tions. As they are exactly solvable their dynamical structure can be inferred from explicit 
expressions. 

Our study of the dynamics of stochastic growth equations on growing domains begins 
with a stochastically forced diffusion equation, known as the Edwards- Wilkinson (EW) equa- 
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tion [10| in this context, which reads 

d t h = DV 2 h + F + £(y,t), (3) 
where £(y, t) is a zero-mean Gaussian white noise which correlation is 

(e(y,W,0) = ^(y-y')^-0, (4) 

D is the diffusion constant, F the constant rate at which mass enters the interface and e the 
noise intensity, all these parameters being positive real numbers. We start considering the 
conservation law in integral form 

j J h(y,t)dy = I [-V-j + ^(y,t)]rfy, (5) 

J St J St 

where St is the uniformly growing domain, j = -DV/j is the current generated by diffusion, 
and JF(y,t) = F + £(y, t) is the EW growth mechanism By applying the Reynolds 

transport theorem we find 



j I h(y,t)dy = I [d t h + V-(vh)]dy, 

J St J St 



(6) 



where v(y, t) denotes the flow velocity generated by the growing domain. Valid as it is for 
any domain, the integral conservation law may be expressed in the local form 

d t h + V- (yh) = DV 2 h + T{y,t). (7) 

In this equation we readily identify two new terms, the advection one v- V/t, and the dilution 
one hV ■ v. For every y G S t , that has evolved from y G S to , we find v(y,t) = dy/dt. Let 
us now concentrate in one- dimensional domains and then move to higher dimensions. In 
this case uniform growth translates into y = g{t)yo, where g(t) is a temporal function such 
that g(t ) = 1. This yields v = yg/g, and thus 

d t h + £ (yd y h + h) = Dd 2 y h + F + £(</, t). (8) 

For a one-dimensional domain (0, L(t)), with L(t) = g(t)Lo, we change the spatial coordinate 
x = yLo/L(t), where Lq = L(to), in order to map the problem into the interval (0, Lq). This 
transformation counterbalances advection, and so the resulting equation reads 



dh fLoY r>^-^ + F + J^( x ,t), (9) 



dt \L(t) J dx 2 g y L{t) 



where we have used the fact that the noise is delta correlated. The dilution term has become 
hV • v = — (g/g)h. Dilution has a transparent physical meaning: as the domain grows the 
incoming mass becomes distributed in a larger (d— dimensional) area. Now we assume that 
the growth function adopts the power law form g(t) = (t/t ) 7 , where the growth index 7 > 0, 
to find 

The growth index 7 is a new degree of freedom of this problem; it cannot be deduced from 
the other model parameters, and has to be measured directly from the physical system under 
study. 

Now we move to a more general situation in which we consider an arbitrary diffusion 
operator of order ( and an arbitrary spatial dimension d, see Eq. ([7]). From now on the 
d— dimensional coordinates will be denoted x — > x and y — > y for simplicity. In this case, 
we can proceed exactly in the same way as in the one-dimensional situation to find, instead 
of Eq. ( Tl0|) . the equation 

o t h = -d ivi^ -Y k+F+ (t)" 772 ^' *)> ( n ) 

where the fractional operator |V|^ is to be understood in terms of the Fourier transform. 
Special values of ( yield some of the well known equations in this topic, as the EW equation 
for ( = 2 [lo| and the Mullins- Herring equation for ( = 4 [12I, [l^. There is still another 
way of deriving an equation similar to (TTTT) but in which the dilution term is not present. If 
we just considered the dilatation transformation x —>■ (t/t^x instead of domain growth we 
would find 

d t h = -D(Jj |V|^ + F+(JJ fat). (12) 
The difference among equations ( TTTT) and ( fl2l) appears already in the amount of mass arriving 



at the interface [14] . In absence of external sources of mass, i. e. F = e = 0, and for no flux 
boundary conditions the total mass on the surface is conserved 

pL(t) pL{t) pL pL 

/ •■• / h(y,t)dy= / ■■• / h(x,t )dx, (13) 
Jo Jo Jo Jo 

for the dilution dynamics (TTTT) . In the no dilution situation corresponding to equation (TT2|) 

we find 

rL(t) rL{t) / ^ \ d-y r-Lo r-L 

/ •■■ / h(y,t)dy=(-) / •■■/ h(x,t )dx. (14) 
^0 ^0 \ r o/ Jo Jo 



This second case, as we have already mentioned, is pure dilatation, which implies that not 
only the space grows, but also the interfacial matter grows at the same rate, in such a way 
that the average density remains constant. Note that this process of matter dilatation, as 
well as the spatial growth, are deterministic processes. 

We have analyzed both types of dynamics, in the absence and presence of dilution, and 
found a number of measurable consequences. Temporal dynamics can be studied by means 
of the temporal auto-correlation 

is (h(x,t)h(x,t')) /min{t,t'}\ A „ . , , . 

where A is the auto-correlation exponent and (-)o denotes the average with the zeroth mode 
contribution suppressed. When dilution is considered the auto-correlation exponent takes 
the form 

A = 

or alternatively 

where /? = 1/2 — d/(2() and 

In absence of dilution we find 

A = < 

or alternatively 
where /? = 1/2 — d/{2() and 

From these formulas one can clearly read that when dilution is suppressed there is no mecha- 
nism for correlations propagation. Otherwise dilution is the responsible for the propagation 
of correlations in the fast growth regime, i. e. when 7 > 1/C,. 



P + d/C if 7 < 1/C, 
P + ld if 7>VC, 



\ + I - d 1 if 7 < 1/C, 
if 7 > 1/C, 



1 — A 

2 2( 



C/(1-7C) if 7<1/C, 
00 if 7 > 1/C- 



(16) 



A = /?+-, (17) 



z x = min{C, I/7}. (18) 



(19) 



X = (3+—, (20) 

ZeS 



(21) 
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Complementary information can be obtained from the interface persistence. The persis- 
tence of a stochastic process denotes its tendency to continue in its current state. When 
considering the dynamics of a fluctuating interface, one refers to the persistence probability 
P + (ti,t 2 ) (P-(ti,t 2 )) as the pointwise probability that the interface remains above (below) 



its profile at t\ up to time t 2 > t% 16j, |17|]. Herein, as in 18j, we concentrate on the case 



in which the initial profile is flat, and we suppress the contribution coming from the zeroth 
mode again. For the stochastic differential equations under consideration the symmetry 



h, 



-h n for all Fourier modes n^O holds, imp 



times t 2 3> t\ we have the power law behavior |16|, Il7 ] 



ying the equality P + = P = P. For long 



P{h,t 2 ) ~ (ti/t 



2) , 



(22) 



defining the persistence exponent 9. It was previously calculated in the limit £ — ► oo when 

7 = o |rj 

„ 1 2y/2-ld 

° « 2 + —2—? ' ' 



up to higher order terms, and in this same limit when d = 1 and 7 = 1 18] 

9 



1 

2 



1 



(24) 



up to higher order terms and in the absence of dilution. We have calculated the persistence 
exponent ll], again in the limit ( — ► 00 and assuming the inequality 7 > l/£, in the 
presence of dilution 

»»5 + *r-|, (25) 



up to higher order terms, and in the absence of it 

9 



1 
2 



d 

2C' 



(26) 



generalizing the previous result [18j . The exponent 9 = 1/2 characterizes neutrally persistent 
interfaces, which are those deprived of a relaxation mechanism (i. e. D = in equation 
(fl2|) ). For 9 < 1/2 the interface is persistent and for # > 1/2 it is antipersistent. Note that 
if dilution acts on the interface then it is antipersistent, as in the case of no domain growth; 
contrarily, if dilution is not present, the interface becomes persistent. 

Before we start calculating spatial correlations let us note that domain growth induces 
the length scale \x — x'\ ~ t( 1 ~C7)/C_ First we show the scaling form that the two points 
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correlation function adopts for "microscopic" spatial scales \x — x'\ <C t^ 1 ^i/C in the fast 
growth regime. In this case one has 14] 



(h(x, t)h(x', t)) - (h(x, t)) 2 « \x - x'\ c - d t^- d) F [\x - x'|t (C7 - 1)/c ] 



(27) 



or in Lagrangian coordinates \y — y'\ = \x — x'lt 7 

(%, t)h(y\ t)) - (h(y, t)) 2 « \y - y'\^ d F 



\y-y 



(28) 



where we have assumed the inequality Q > d and the statistical isotropy and homogeneity of 
the system in the limit in which the scaling form holds. As dilution does not act on such a 
microscopic scale, these results are independent of whether we contemplate dilution or not. 
Things are different for macroscopic length scales \x — x'\ ^> t^ 1 ^^^. In this limit dilution 
has a measurable action, and when it is included in the interface equation of motion the 
resulting correlation is llj 



(h(y,t)h(y',t))-(h(y,t)} 2 ~t6(y-y') 



(29) 



which is simply the short time limit of equation (|28jl . If dilution is suppressed we find 



however 



20 



21 



(h(y,t)h(y',t)) ~ { 



tS(y-y') if 7<1M 

t\n(t)5(y-y') if >y = l/d, 
V d 5(y-y') if 7 > l/d. 



(30) 



In this case we see that 
dependent prefactor 



20 



or 



ast enough growth memory effects appear and modify the time 
211 ] . The increase of this prefactor reflects the mass excess that 
enters the interface when dilution is not operating as shown in equation f[T4l) . A consequence 
of all these correlations is the scale dependent fractal dimension 



df(\x — x'\,t) 



1 + (3d -C)/2 if 
d+ 1 if 



\x - x'\ < t^-CTVC, 
\x - x'\ > t^-frrVC 



(31) 



which is independent of whether we contemplate dilution or not, and reveals the interface 
multifractality. These asymptotic values suggest the self-similar form of the fractal dimension 



df = d 



f 



\x — x 



£(1-C7)/C J ' 



(32) 
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which would imply its invariance with respect to the dilatation x — > bx, t — > b Zf t, and 
df — > b a fdf, for 2/ = £/(l — £7), a/ = and 6 a real number strictly greater than one. 

As a final note let us mention that the assumption ( > d is fundamental in order to get 
the correlations specified by equation (1251) . As we have seen, in this case the dynamical 
exponent is universal and given by z = (. For 7 = and ( = d = 1 the dynamical exponent 
is still universal and given by z = 1; however for 7 = £ = d = 1 and in the absence of 



dilution this exponent becomes non-universal and given by z = F/(D + F) G (0,1) 19] 



The presence of dilution restores universality and the non-growing domain result z 



l4. 



III. NONLINEAR GROWTH 

The open questions in this topic are related, not surprisingly, to the appearance of non- 
linear terms in the corresponding equations of motion. One of the most popular nonlinear 
models in this context is the KPZ equation, which as we have commented in the Introduc- 
tion is related to the biologically motivated Eden model. As we will see, understanding the 
KPZ equation on a growing domain may shed some light on some of the properties of the 
classical version of this model. 

The KPZ equation on a growing domain reads [ijj] 

dth = u (£) ^ V 2 /> + \ (*) ^ iyhf -^h + jFt^ + (I) fr, t). (33) 
Of course, if we just considered the dilatation x — > {t/t^x we would find 

dth = v (I) ^ V*h + \ (I) ^ {Vhf + 7 FT-i + (|) ^ £(*, t). (34) 

As we have shown in the previous section, the dilution mechanism fixes the Family- Vicsek 
ansatz in the fast growth regime. In the radial Eden model case, assuming it belongs to 
the KPZ universality class, we would have z — 3/2 in d = 1 and 7 = 1. And so, one 
would nai'fly expect that the resulting interface is uncorrelated and we have to resort on 
dilution effects in order to fix the Family- Vicsek ansatz and get rid of memory effects. 
But here comes the paradoxical situation. There are two main symmetries associated with 
the (i-dimensional KPZ equation: the Hopf-Cole transformation which maps it onto the 



noisy diffusion equation 22| and Galilean invariance which have been traditionally related 



to the non-renormalization of the KPZ vertex at an arbitrary order in the perturbation 



expansion 23||. In the case of the no-dilution KPZ equation (I34p both symmetries are 
still present. Indeed, this equation transforms under the Hopf-Cole transformation u = 
exp[Xh/(2u)} to 



d t u = v\ — \ V u + -^—t u + — 



ft \ dj/2 



(35) 



which is again a noisy diffusion equation and it can be explicitly solved in the deterministic 
limit e = 0. We find in this case 



u(x, t) 



(1 -2 7 ) d / 2 exp[FAr/(2z/)] 



which corresponds to 



exp 



h(x, t) 



2v J (1 - 2 7 ) d / 2 exp[FAt 7 /(2z/)] 



\x-y\\l-2 1 ) 
|x-y| 2 (l-2 7 ) 



u(y,t )dy, (36) 



A 



h(y,to) 



dy 
(37) 

for given initial conditions u(x, to) and /i(x, to)- It is clear by regarding this formula that 
decorrelation at the deterministic level will happen for 7 > 1/2. It is still necessary to 
find out if at the stochastic level this threshold will be moved to 7 > 2/3. If we consider 
the dilution KPZ equation (13"3l then transforming Hopf-Cole we would find the nonlinear 
equation 



,'t \ 27 2 d l i / \ 1 7^.7-1 , A {to 



d-y/2 



£{x,t)u, 



(38) 



which may be thought of as a time dependent and spatially distributed version of the Gom- 
pertz differential equation 24j . In this case it is not evident how to find an explicit solution 
at the deterministic level and what would be its decorrelation threshold. 
Galilean invariance means that the transformation 

A 



x 



x — Aft, 



h — ► h + vx, 



F v 



(39) 



where v is an arbitrary constant vector field, leaves the KPZ equation invariant. In case of 
no dilution this transformation can be replaced by 



A 



x — ► x 



h + vx, F -> F - —v 2 t 2 ^t^ 3 \ (40) 

27 



1 - 2 7 

which leaves invariant equation (|34l) . If we consider dilution, then it is not clear how to 
extend this transformation to leave equation fl33l invariant. The main difficulty comes 
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from the dilution term which yields a non-homogeneous contribution to the dynamics as a 
response to the transformation h —>■ h + vx. So in summary we may talk of a certain sort of 
Galilean invariance which is obeyed by the no-dilution KPZ dynamics (1341) and is lost when 
dilution is taken into account. If it were found that the dilution equation (J33D obeys the 
traditional KPZ scaling (at least in some suitable limit), then that would put into question 
the role that Galilean invariance has in fixing the exponents. The KPZ critical exponents 
are believed to obey the scaling relation a + z = 2 in all spatial dimensions, a relation that 



has been traditionally attributed to Gali 



been recently put into question 



26 



27 



ean invariance, although this interpretation has 



There is still another fundamental symmetry of the KPZ equation, but this time it just 
manifests itself in one spatial dimension: the so called fluctuation-dissipation theorem [5, 

n 

16J. It basically says that for long times, when saturation has already being achieved, the 
nonlinearity ceases to be operative and the resulting interface profile would be statistically 
indistinguishable from that created by the EW equation. For fast domain growth, we know 
from the linear theory that the interface never becomes correlated, and it operates, in this 
sense, as if it were effectively in the short time regime a consequence, 

the fluctuation-dissipation theorem is not expected to play any role in this case. Of course, 
this result would be independent of whether we contemplated dilution or not. 

IV. DISORDER 

It is quite natural to consider propagation into disordered media in the context of biolog- 
ical growth. One could, for instance, imagine the development of a bacterial infection inside 
a host body. This is, of course, the propagation of a bacterial front inside a medium with a 
extremely low degree of symmetry. 

Apart from the classical KPZ equation which is driven by thermal noise, a different 
version in which this noise is replaced by quenched disorder has been considered in the 
literature 0, ^ 

d t h = vV 2 h + ^(Vh) 2 + F + ^(x,h), (41) 
(£{x,h)) = 0, (42) 
(£(ar, h)£(x', h')) = eS(x - x')5{h - h'), (43) 
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which is nonlinear even when A = 0. All the problems considered in this section are open in 
this simpler case as well. A way of understanding this equation is considering the simplified 
random deposition version of it 

d t h = F + £(x,h), (44) 

which is actually an ordinary differential equation in which the position x acts just as a 
parameter, at least if the quenched disorder is conveniently regularized. And thus, let us 
consider the auxiliary problem 



dh 
~dt 

(rj(h)) = 0, (46) 



F + r)(h), (45) 



(r]{h)r]{ti)) = eS(h-h'). (47) 
It is similar to the stochastic problem 



f = F + V {t), (48) 
(V(t)) = 0, (49) 
( V (t)v(t')) = e5(t-t'), (50) 



which solution is 

h{t) = h(t ) + F(t - t ) + V~t W(t - t ), (51) 

where to is the initial time and W(t) is a Wiener process, so we basically have two super- 
posed motions: constant drift and Brownian motion. Classical Brownian motion describes 
a particle choosing its direction of motion randomly every time step. If the noise is position 
dependent, i. e. r\ — rj(h), then the direction of motion is already prescribed in every spatial 
point. For a discrete version of this process, say an unbiased random walk on Z for which 
the jump direction is specified in every site, if the system revisits any location, then it is 
trapped forever (jumping forth and back in the last two visited sites). The only way to 
prevent (this of other sort of) trapping is to consider a sufficiently large F, so the system 
is evolving over new positions all of the time. Something similar happens in equation fj4Tl) : 
for values of F smaller than a critical one F r the interface becomes pinned, while for larger 

n 

values the interface propagates [5|]. 

It would be interesting to analyze the interplay of spatial disorder with a growing domain 
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size. We again have two possibilities, the equation with dilution 

d t h = v (|) ^ V 2 /. + \ ^ (Vhf -^h + jFt^ + (f) ^ /i), (52) 

and the one without it 

dth = u (^j ^ V 2 /. + ^ ^ (V/,) 2 + 7-Ff y ~ 1 + ^ fc). (53) 

Note that, contrary to the classical case in which the interface propagates linearly in time, 
in this case the velocity of propagation would be ~ P '. It would be interesting to clarify 
whether 7 = 1 plays a critical role in the dynamics or not. In other words, whether the 
interface is always pinned for 7 < 1 and always moving for 7 > 1, or if there is a dependence 
on the parameters values in these cases too. If the second situation held, it would be yet 
necessary to clarify whether there are other possible critical values for 7. According to the 
microscopic description commented in this section, it seems plausible that 7 = 1 is indeed 
the critical value of the growth index separating pinned and unpinned regimes. 

Even when 7 = 1 there is an interesting question associated with dilution. We know 
that dilution keeps constant the amount of matter on the interface, while suppressing it we 
get a mass excess. For linearly in time growing interfaces we know that the value of F, 
which describes the amount of matter arriving at the interface, controls the possibility of 
interface pinning/unpinning. If no dilution is present, the mass excess could act as effectively 
increasing the value of F, and thus facilitating interface unpinning. It would be interesting 
to quantitatively determine how much the threshold of pinning is moved in the absence of 
dilution, if this is indeed the case. 



V. SUMMARY AND CONCLUSIONS 



As we have seen, a consequence of the linear theory is that dilution erases the memory 



effects and this way restores the classical Family- Vicsek ansatz |20l. |21|. Otherwise, for fast 
domain growth, a series of unexpected consequences arise, as the modification of the random 
deposition correlation, the lost of the antipersistent character of the fluctuating interface 
and even the appearance of non- universal critical exponents. In this respect, dilution can 
be thought of as the mechanism which maintains some of the most characteristic features of 
surface growth when we let the domain size grow in time. 
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On the other hand, some of the well known symmetries of the KPZ equation, as Galilean 
invariance and mapping to the directed polymer problem, are maintained in a dilating set- 
ting but lost by virtue of dilution. So, in principle, one would expect that in the absence of 
dilution memory effects could be present in the KPZ dynamics, and this way some its char- 
acteristic features would be lost. However, it is somehow paradoxical that it is exactly this 
absence of dilution what maintains the classical symmetries of this equation. Furthermore, 
these symmetries have been sometimes considered as necessary ingredients in the resulting 
KPZ scaling. If we found that the KPZ equation in a rapidly growing domain and in pres- 
ence of dilution behaved in a similar way to its classical counterpart, that would suggest 
that the symmetries present in the standard situation are not playing such a necessary role. 

As we have already mentioned, the motivation for studying radial growth models such 
as the Eden or different ones partially comes from the possible similarity of these with 
some forms of biological development, such as for instance bacterial colonies formation. The 
results of our study can be translated into this context to obtain some simple conclusions, 
provided the modelling assumptions make sense for some biological system. The structure 
of a rapidly developing bacterial colony would be dominated by dilution effects, originated 
in the birth of new cells which volume causes the displacement of the existent cells. If the 
rate of growth is large enough this motion will dominate over any possible random dispersal 
of the bacteria. It is remarkable that such a consequence simply appears by considering 
domain growth, while it is not necessary to introduce corrections coming from the finite size 
of the constituents. This is the dilution dominated situation we have formalized by means 
of the (decorrelation) inequality 7 > If we were to introduce some control protocol 

in order to keep the consequences of bacterial propagation to a minimum we would need 
to eliminate colony constituents (possibly randomly selected) at a high enough rate so the 
effective growth velocity were one that reversed the decorrelation inequality. For the two 
dimensional radial Eden model, accepting it belongs to the KPZ universality class, one finds 
7 = 1 and z = 3/2. If z played the same role for the nonlinear KPZ equation as ( for the 
linear equations considered herein (as it is reasonable to expect), the Eden model would 
be in the dilution dominated regime. In order to control it we would need to eliminate its 
cells at rate such that the effective growth rate obeyed 7 < 2/3. For the three dimensional 
Eden model, if its behavior were still analogous to that of the KPZ equation, we would find 
z > 3/2 and thus a greater difficulty for control. Note that for the particular growth rules of 
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the Eden model one would need to eliminate peripheral cells in order to control the system. 
This would not be so in the case of an actual bacterial colony, for which bulk cells are still 
able to reproduce, and so cell elimination could be performed randomly across the whole 
colony Of course, these conclusions are speculative as long as stochastic growth equations 
are not proved to reasonably model some biological system. With respect to the problem of 
the experimental verification of the theoretical results, one would be interested in finding a 
method for measuring the exponents. These could perhaps be measured using the long time 
dependence of the interface variance on the initial system size and time 



that we found for both the dilution and no-dilution linear dynamics in [11]. This way one 
could in principle experimentally determine both exponents a and /3 and as a consequence 
the ratio z = a/ (3 in the fast growth regime. Subtracting the inverse of the so obtained value 
for the dynamic exponent z from the measured value of the growth index 7, which should be 
easily obtainable from experiments, one could estimate the distance from correlation, and 
in turn the possible necessary strength of a control protocol. 
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